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Abstract 

£vq ■ The effects on the 4/ and 5d levels of Ce 3+ of its first coordination shell geometry in Ce-doped 

oxides with a Z?2 8-fold site, like garnets, are studied with embedded cluster, wave function based ah 

O 

initio methods. The only deformations of a CeOs cube that are found to shift the lowest 4/ — > 5d 

transition to the red are the symmetric Ce-0 bond compression and the tetragonal symmetric 

bond bending. These results are analyzed in terms of centroid and ligand field stabilization energy 

c3. 

r^ differences. The splittings of the upper 5d levels and of the 4/ levels are also discussed. 
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I. INTRODUCTION 

Yttrium aluminum garnet Y3AI5O12 (YAG) doped with Ce 3+ is a well known phosphor 
with a blue Ce 3+ 4/ — y 5d absorption and a corresponding yellow 5d —y 4/ emisssion, 1 which 
is used in InGaN blue LED based, white light solid-state lighting devices. 2 ' 3 In the search for 
alternative phosphors with the best efficiencies and tailored color-rendering indexes, doping 
Ce 3+ in other artificial garnets has been a line of work. It led, for instance, to the discoveries 
of the Lu2CaMg 2 Si30i2:Ce 3+ orange phosphor 4 and the Ca3Sc2Si30i2:Ce 3+ green phosphor. 5 

Most garnets can be described in terms of a 160 atom body-centered cubic unit cell (80 
atom primitive cell) of the Ia3d (230) space group, which contains eight formula units of 
A3B2B3O12. A, B' and B" are cations of different nominal charges in different symmetry 
sites. In Y3AI5O12, for instance, A=Y, B'=A1, and B"=A1, 6 all with oxidation states +3. 
In Ca3Sc2Si30i2, A=Ca, B'=Sc, and B"=Si, 7 with respective oxidation sates +2, +3, and 
+4. In Lu2CaMg2Si30i2, 2/3 of the A sites are occupied by Lu and 1/3 by Ca atoms, 4 and 
B'=Mg, and B"=Si, with respective oxidation states +3, +2, +2 and +4. In all cases, A 
occupy 24(c) sites of 8-fold coordination, B' 16(a) sites of 6- fold coordination, and B" 24(d) 
sites of 4-fold coordination, all of them at fixed positions, with the remaining 96 O atoms 
in (h) sites, which depend on three x, y and z internal parameters. 6 Optically active Ce 3+ 
impurities substitute for A atoms at (c) sites, which have the local symmetry of the D2 point 
group and a first coordination shell made of 8 oxygen atoms. 

The energies of the 4/ and 5d levels of Ce 3+ in particular garnets and in other oxides 
depend on the bonding and electrostatic interactions between Ce and the hosts, and they 
can be calculated in a one-by-one basis, with reasonable accuracy, by means of ab initio 
methods, both in D2 local symmetry like in Ce-doped YAG 8 and in lower local symmetries 
like in Ce, La-doped YAG 9 and Ce,Ga-doped YAG. 10 These energies are dominated, up 
to first-order, by the bonding interactions between Ce and its first oxygen coordination 
shell, subject to the basic confinement embedding effects of the host. Other specific host 
embedding effects are important at higher orders of approximation and they refine the 
results. 

This paper is aimed at providing a basic study of the energies of the levels of the 4/ 1 
and bd 1 manifolds of Ce 3+ in D2 8-fold oxygen coordination. We report the results of an 
ab initio calculation of them after consideration of the bonding interactions with the first 



oxygen coordination shell only, under the effects of a cubic, non host specific, confinement 
embedding potential. The present results are a common reference for these levels of Ce 3+ in 
specific garnets (and other oxides with 8-fold coordination), since the latter can be considered 
as resulting from the former after including the host specific embedding effects. 

II. METHOD 

In order to know the 4/ 1 and 5d l energy levels of Ce 3+ under the effects of its interac- 
tions with a coordination shell of eight oxygens in a D 2 symmetry site of a solid oxide, we 
calculated the corresponding energy levels of the (CeOs) 13- cluster embedded in a Oh host 
potential. We will describe later the details of the ab initio wave function based quantum 
mechanical calculation. With he choice of a cubic host embedding potential, all the energy 
splittings and changes due to D 2 distortions are ascribable to the interactions between Ce 
and the first-shell O atoms. 

It is convenient to look at the D 2 energy levels as derived from the more familiar Oh 
levels after geometry distortions and this is what we did in this work. In consequence with 
this, we decided to take a cubic CeOs Oh moiety as a reference and to chose six atomic 
displacement coordinates Si — Sq that transform according to irreducible representations of 
the Oh point symmetry group for the six degrees of freedom of the CeOg D 2 moiety. They 
are defined in Table I and represented in Fig. 1. Si represents the breathing (or symmetric 
bond stretching) of the reference cube and it transforms according to the totally symmetric 
irreducible representation oi ff ; S 2 represents an asymmetric bond stretching of the crosses 
ai — cl 2 — 03 — 04 and ei — e% — e 2 — e 4 and it transforms according to the e g e sub-especies of the 
doubly degenerate e g irreducible representation; 11 S 3 and S 4 are the e g 9 and e g e symmetric 
and asymmetric bond bendings of the a and e crosses; and S5 and Sq are the e u and e u e 
symmetric and asymmetric twistings of the crosses. 

The details of the quantum mechanical calculation are the following. We performed 
ab initio calculations on a (CeOs) 13- cluster that include Ce-0 bonding effects, static 
and dynamic electron correlation effects, scalar and spin-orbit coupling relativistic effects, 
and cubic host embedding effects. For each D 2 nuclear configuration of the (CeOs) 13- 
embedded cluster, we performed, in a first step, scalar relativistic calculations with the 
many-electron second-order Douglas-Kroll-Hess (DKH) Hamiltonian. 12 ' 13 These were state- 



average complete-active-space self-consistent-field 14-16 (SA-CASSCF) calculations with the 
active space that results from distributing the open-shell electron in 13 active molecular 
orbitals with main character Ce 4/, 5d, 6s, which provided occupied and empty molecular 
orbitals to feed subsequent multi-state second-order perturbation theory calculations (MS- 
CASPT2), 17 ~ 20 where the dynamic correlation of 73 electrons (the 5s, 5p, 4/, 5d, 6s electrons 
of Ce and 2s, 2p electrons of the O atoms) were taken into account. In a second step, 
spin-orbit coupling effects were included by adding the AMFI approximation of the DKH 
spin-orbit coupling operator 21 to the Hamiltonian and performing restricted-active-space 
state-interaction spin-orbit calculations (RASSI-SO) 22 with the previously computed SA- 
CASSCF wave functions and MS-CASPT2 energies. 

In all these calculations, the Hamiltonian of the (CeOg) 13- cluster is supplemented with 
a cubic embedding Hamiltonian. For this we chose the ab initio model potential embedding 
Hamiltonian (AIMP) 23 of a SrO host with CsCl lattice structure (Pm3m, no. 221), which 
provides a perfect 8-fold cubic oxide coordination of the cationic site. We used a lattice 
constant a = 2.87 A (the one that gives a Ce-0 equilibrium distance in the cubic (CeOs) 13- 
embedded cluster of 2.34 A at the MS-CASPT2 level of calculation, which is the Ce-0 
distance of the reference CeO§ cube in an hypothetical undistorted Ce-doped YAG 6 ). The 
embedding potential of this 8-fold coordinated cubic oxide was computed in Hartree-Fock 
self-consistent embedded- ions calculations (SCEI), 24 in which the input embedding AIMPs 
used for the Sr 2+ and O 2- ions of SrO are consistent with the output AIMPs obtained out 
of the HF orbitals of the embedded Sr 2+ and 2 ~ ions. 

III. RESULTS 

The dependence of the 4/ — > 4/ and 4/ — > 5d transitions on the D 2 distortions of the 
first coordination shell are shown in Fig. 2 without spin-orbit coupling and in Fig. 3 with 
spin-orbit coupling interactions included. The chosen ranges of the Si — 5*6 coordinates span 
their usual values in natural and artificial garnets, which are shown in Table II and in the 
figures. Besides the transitions, the figures include two additional lines: One is the difference 
between the bd and 4/ energy centroids (or average energies of the five bd 1 levels and the 
seven 4/ 1 levels), AE centroidJd = E centroid (hd l ) - E centroid (Af l )\ the other is the difference 
between the ligand field stabilization energies of the hd and 4/ lowest levels AEu gand f ie i d)d f = 



[E centroid (hd l ) - E^d 1 )} - [E centroid (4f l ) - E^Af 1 )}. The fact that the lowest 4/ -»• hd 
transition equals their subtraction, Eiihd 1 ) - Ei(4f l ) = AE centroid j d - AE ligand f ield4 f, has 
been used to analyze the reasons behind red and blue shifts of this transition. 9,10 

Let us first focus on the lowest 4/ — > 5d transition, whose reverse is the luminescence of 
Ce 3+ in garnets and other oxides. It is clear that only two coordinates have an important 
impact on it: the breathing mode Si, which does not distort the cubic symmetry, and the 
symmetric bond bending mode 5*3, which gives a e g 9 tetragonal distortion of the CeOs cube. 

Symmetric compression of the Ce-0 bonds lowers the 4/ — ¥ 5d transition. This is almost 
entirely due to the increase of the 5d ligand field it produces, which lowers the first 5d l level 
with respect to its centroid. This effect is slightly compensated with a small increase due 
to the centroids: bond compression slightly shifts the bd centroid upwards with respect to 
the 4/ centroid. The latter observation has been made before 10 ' 25 and it contradicts the 
predictions of the usual model for the centroid energy difference, 26 ~ 28 according to which 
it should be smaller for smaller bond lengths; nevertheless, the model has been useful to 
rationalize 5d — > 4/ luminescences. 29 The contribution of the breathing mode S\ to the 
lowest 4/ —7- 5d transition represented in Fig. 3 fits 8255 S\ — 4803 S 2 ; for cubic CeOs 
moieties Si = 2y2 Ad, Ad being the Ce-0 bond length change with respect to the reference 
cube; in our case the reference cube has Ce-0 bond length of 2.34 A and a lowest 4/ — > 5d 
transition of 26000 cm -1 . 

Symmetric bond bending along the S3 coordinate lowers de 4/ — » 5d transition. Among 
the D2 distortions of the cube, it is the only tetragonal one. Also, it is the only one with 
a significant impact on the transition (for the ranges of S2 — Sq shown by garnets). The 
enhancement of tetragonal distortions of the cube as responsible for the red-shifts of the 
Ce 3+ luminescence has been suggested by Cheetham and coworkers 30 and it is supported by 
the present investigation. The image emerging from an analysis in terms of configurational 
energy centroids and ligand field stabilization is one in which the lowering of the first 5c? 1 
level with respect to the ground state is due in almost equal amounts to an increase of 
the ligand field (of its tetragonal component) and to a reduction of the centroid energy 
difference. Again, the latter effect contradicts the prediction of the Judd-Morrison model: a 
very small increase of the centroid energy difference out of a very small increase of the Ce-0 
distances. The contribution of the tetragonal mode S3 to the lowest 4/ — > 5d transition 
represented in Fig. 3 fits 1844 S 3 - 2885 S| for S 3 < and -1876 S 3 - 141 Si for S 3 > 0. 



This contribution is due in almost equal amounts to the lowering of the centroid energy 
difference and the increasing of the ligand field. The centroid effect fits 10 S3 — 2695 S 3 for 
S 3 < and 12 S 3 - 2581 Sf for S 3 > 0; the ligand field effect fits 1834 S 3 - 190 S| for S 3 < 
and -1887 S 3 + 2440 Si for S 3 > 0. 

Apart from Si and S3, only the asymmetric stretching S2 has an effect, although very 
small, on the first 4/ — > 5d transition (it fits —2033 Sf in the range of the figures). 

Regarding higher 4/ — > 5d transitions, several remarks can be made from Figs. 2 and 
3. The gap between the first and the second bd levels related with the cubic 2 E g is due 
to the symmetric bending S3 and to the symmetric twisting S 5 , with a significantly higher 
contribution from the latter in garnets; however, different splittings in different garnets 
would be due to S3. The splitting of the three upper levels, mostly related with the cubic 
2 T2 9 , comes basically from the symmetric bending S3 and the asymmetric stretching S2, 
with a small contribution from the symmetric twisting S5 and negligible contributions from 
the asymmetric bending and twisting S4 and Sq. The larger Si and the absolute values of 
S2 and S3, the higher probability for the third 5d level to appear below the conduction band 
of the garnet. 

Finally, the 4/ levels are shown with more detail in Figs. 4 and 5, without and with 
spin-orbit coupling respectively. These may be interesting because the 5d — > 4f emission 
is made of the superposition of the seven individual emissions and the full with and shape 
of the emission band depends, although not only, on the relative positions of the 4/ levels. 
Since the emission to the highest 4/ level has a minor contribution to the full emission band 
shape, 8 we can have a predictive hint by looking at the six lowest 4/ levels. According 
to Fig. 5, the width seems to be controlled by S2 and S3, so that the garnets with higher 
absolute values of these two coordinates would tend to have wider emission bands because of 
their higher 4/ level separation. Of course, the width of each individual 5d — > 4/ emission 
will depend on the 5d and 4/ equilibrium offsets. 

IV. CONCLUSIONS 

The levels of the 4/ 1 and bd 1 manifolds of Ce 3+ in D2 8-fold oxygen coordination have 
been calculated ab initio as a function of D2 deformation coordinates of a reference cube, 
with the goal of pinpointing the effects of the geometry of the first coordination shell. A 
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(CeOs) 13 cluster under the effects of a confinement AIMP embedding potential of cubic 
symmetry has been used at a level of calculation all-electron DKH for the Hamiltonian and 
SA-CASSCF/MS-CASTP2/RASSI-SO for the wave functions. The results include bonding, 
correlation, scalar relativistic, spin-orbit coupling, and embedding interactions. 

It is found that the lowest 4/ — > 5d transition shifts significantly to the red only as a 
consequence of symmetric Ce-0 bond compression and tetragonal symmetric bond bending. 

The increase of the 5d cubic ligand field dominates the effects of the bond compression, 
whereas the effect of the tetragonal bond bending distortion is divided in almost equal 
amounts between an increase of the 5d tetragonal ligand field and a reduction of the energy 
difference between the 5d and 4/ energy centroids. 

The splittings of the upper 5d levels and of the 4/ levels have also been studied. 
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TABLE I: Definitions of the D^ totally symmetric displacements of the CeOs moiety chosen in 
this work: Stretching, bending, and twisting Oh symmetry coordinates. The labels of the oxygen 
atoms and the chosen cartesian axes are defined in Fig. 1: In the reference cube, the symmetry 
independent oxygen atoms a\ and e\ are located at (v2, 0, l)d/v3 and (0, \/2, l)d/y/?> respectively, 
with d being the Ce-0 distance. R is the Z?2 symmetrization operator: the normalized addition 
of the group symmetry operations (the identity and the three 180° rotations around the cartesian 
axes). 8x a \ is the x cartesian displacement of the Oxygen atom al from its position in the reference 
cube: 5x a \ = x a i — x a \ re t\ identical definitions stand for the other cartesian displacements. 





R = 


2 ( I + @2x + C<ly + C2z) 


Displacement 


Oh irrep 


Definition 


symmetric stretching 


°is 


Sl = 7e R 


yV^Sx a i + 5z a ij + [y2Sy e i + 5z el j 




asymmetric stretching 


e g e 


S2 = 7e R 


yV^Sx a i + 5z al J - (y26y e i + 5z el J 




symmetric bending 


e g 9 


53 = 75^ 


(-5x a i + V2 5z a ij + [-5y e i + \/2 8z e ij 


asymmetric bending 


e g e 


Si = 7s R 


{-5xai + V2 5z a ij - y-Sy el + y/2 5z el J 


symmetric twisting 


e u 9 


S$ = ^ R [Sy a i + Sx el ] 


asymmetric twisting 


e«e 


Se = 73 R ^ Val + Sxe1 ^ 
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TABLE II: Crystallographic data of some garnets: lattice constant a and special position (d) of the Ia3d (230) space group x, y, z. Local 



Z?2 8-fold oxygen coordination of the fixed position (c) occupied by Ce as a dopant: cation-oxygen distance in the reference cube d re f and 
-D2 displacements Si - Sq as defined in Table I; Si is defined with respect to a cube with cation-oxygen distance d =2.34 A (note that the 
use of d re f, Si = 0, and 52 — 5*6 or d re f =2.34 A, and Si — S$ are two valid alternatives to define the oxygen coordinations), a, d re f and Si 
- Sq are given in A and x,y,z in fractional cell units. The lowest 4/ — > 5d transitions calculated for the (CeOs) 13- clusters, embedded in 
a common cubic confinement potential and having the local structures of the undistorted garnets, are given in the last column in cm" 1 ; its 
stabilization with respect to a reference cube of d =2.34 A is given in parenthesis. 



Garnet 



Ref. 



.c 



y 



z 



d 



ref 



Si 



So 



s 3 



Si 



S5 Sq 



4/^5d 



Pyrope Mg3Al 2 Si 3 0i2 31 11.4566 -0.03290 0.05010 0.15280 2.243 -0.2731 -0.2191 -0.0768 -0.1215 0.9900 -0.0892 22750 

Almandine Fe 3 Al 2 Si30i2 32 11.5250 -0.03401 0.04901 0.15278 2.267 -0.2070 -0.2334 -0.1126 -0.1279 1.0084 -0.1028 23380 

Spessartine Mn 3 Al 2 Si 3 0i2 33 11.6190 -0.03491 0.04791 0.15250 2.295 -0.1278 -0.2490 -0.1414 -0.1337 1.0281 -0.1243 24100 

Grossular Ca 3 Al 2 Si 3 0i 2 31 11.8515 -0.03760 0.04540 0.15130 2.368 0.0790 -0.2752 -0.2203 -0.1378 1.0725 -0.1909 25670 

LuAG Lu 3 Al 5 0i2 6 11.9060 -0.02940 0.05370 0.15090 2.303 -0.1042 -0.1683 0.0555 -0.0866 0.9766 -0.1044 24910 

YbAG Yb 3 Al 5 0i2 6 11.9310 -0.02960 0.05290 0.15040 2.313 -0.0759 -0.1819 0.0497 -0.0906 0.9858 -0.1286 25140 

Lu 2 CaMg 2 Si 3 0i2 4 11.9750 -0.03510 0.05380 0.15780 2.333 -0.0200 -0.1522 -0.1619 -0.1139 1.0153 0.0956 25240 

ErAG Er 3 Al 5 0i2 34 11.9928 -0.03039 0.05124 0.14915 2.339 -0.0039 -0.2047 0.0297 -0.0944 1.0046 -0.1854 25710 

YAG Y 3 Al 5 0i2 6 12.0000 -0.03060 0.05120 0.15000 2.339 -0.0019 -0.2094 0.0156 -0.1018 1.0099 -0.1613 25680 

Andradite Ca 3 Fe 2 Si 3 0i 2 7 12.0578 -0.03940 0.04870 0.15540 2.398 0.1648 -0.2102 -0.2895 -0.1260 1.0718 -0.0351 26210 

GdAG Gd 3 Al 5 0i2 6 12.1130 -0.03110 0.05090 0.14900 2.368 0.0788 -0.2060 0.0101 -0.0932 1.0180 -0.1957 26340 

LuGG Lu 3 Ga 5 0i2 6 12.1880 -0.02520 0.05700 0.15060 2.322 -0.0520 -0.1378 0.1961 -0.0727 0.9580 -0.0755 25870 

YbGG Yb 3 Ga 5 0i2 6 12.2040 -0.02590 0.05630 0.15190 2.328 -0.0336 -0.1561 0.1596 -0.0899 0.9744 -0.0459 25860 

Ca 3 Sc 2 Si 3 0i2 7 12.2500 -0.04004 0.05010 0.15887 2.427 0.2454 -0.1972 -0.3362 -0.1401 1.0894 0.0842 26580 

YGG Y 3 Ga 5 0i2 35 12.2730 -0.02740 0.05460 0.14930 2.363 0.0647 -0.1629 0.1338 -0.0757 0.9876 -0.1440 26640 

HoGG Ho 3 Ga 5 Oi2 36 12.2900 -0.02740 0.05520 0.15020 2.362 0.0616 -0.1550 0.1293 -0.0771 0.9861 -0.1101 26610 

DyGG Dy 3 Ga 5 0i2 36 12.3060 -0.02780 0.05490 0.14990 2.369 0.0823 -0.1558 0.1194 -0.0752 0.9896 -0.1229 26740 

TbGG Tb 3 Ga 5 0i2 37 12.3474 -0.02820 0.05470 0.14950 2.381 0.1166 -0.1538 0.1115 -0.0712 0.9933 -0.1377 26980 

GdGG Gd 3 Ga 5 0i2 37 12.3829 -0.02890 0.05420 0.14940 2.394 0.1538 -0.1578 0.0901 -0.0719 1.0019 -0.1473 27180 

SmGG Sm 3 Ga 5 0i2 37 12.4361 -0.02920 0.05300 0.14900 2.412 0.2031 -0.1816 0.0777 -0.0818 1.0190 -0.1748 27410 



3250) 
-2620) 
■1900) 
-330) 
1090) 
-860) 
-760) 
-290) 
-320) 
+210) 
+340) 
-130) 
-140) 
+580) 
+640) 
+610) 
+740) 
+980) 
+1180) 
+1410) 



Figure captions 

FIG. 1: D% totally symmetric displacements of the CeOs moiety: Stretching, bending, and twisting 
Oh symmetry coordinates S\ {a\ g symmetric stretching -breathing-), 52 (e g e asymmetric stretch- 
ing), S3 (e g 9 symmetric bending), S4 (e g e asymmetric bending), S5 (e u 9 symmetric twisting), and 
Sq (e u e asymmetric twisting). See Table I for the detailed definitions. 

FIG. 2: 4/ (dashed lines) and 5d (full lines) energy levels of the (CeOs) 13- embedded cluster 
(relative to the ground state) as a function of the S\ - Sq D 2 oxygen displacement coordinates, 
calculated without spin-orbit coupling. The differences between the 5d and 4/ energy centroids 
(dot-dashed lines) and between the ligand field stabilization energies of the 5d and 4/ lowest levels 
(doted lines) are also shown; the lowest 4/ — > 5d transition equals the subtraction of these two. 
Experimental values of S\ - Sq of 21 pure garnets are shown as small vertical lines. 

FIG. 3: 4/ and 5d energy levels of the (CeOs) 13- embedded cluster as a function of Si - Sq, 
calculated with spin-orbit coupling. All levels are D' 2 T§ Kramer doublets. See Fig. 2 caption. 

FIG. 4: 4/ energy levels of the (CeOs) 13- embedded cluster as a function of S\ - Sq, calculated 
without spin-orbit coupling. Dotted lines: 2 A levels; full lines: 2 B\ levels; dashed lines: 2 B 2 and 
2 i?3 levels. Experimental values of S\ - Sq of 21 pure garnets are shown as small vertical lines. 

FIG. 5: 4/ energy levels of the (CeOs) 13- embedded cluster as a function of S± - Sq, calculated 
with spin-orbit coupling. All levels are D' 2 T§ Kramer doublets. Experimental values of S± - Sq of 
21 pure garnets are shown as small vertical lines. 
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Figure 1. Seijo and Barandiaran 
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Figure 2. Seijo and Barandiaran 
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Figure 3. Seijo and Barandiaran 
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Figure 4. Seijo and Barandiaran 
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Figure 5. Seijo and Barandiaran 



